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Abstract 



We obtain a black hole solution in the Einstein-Gauss-Bonnet theory for the string cloud model in a 
five dimensional spacetime. We analyze the event horizons and naked singularities. Later, we compute 
the Hawking temperature Tk, the specific heat C, the entropy S, and the Helmholtz free energy F of the 
black hole. The entropy was computed using the Wald formulation. In addition, the quantum correction 
to the Wald's entropy is considered for the string cloud source. 

We mainly explore the thermodynamical global and local stability of the system with vanishing or 
non-vanishing cosmological constant. The global thermodynamic phase structure indicates that the 
Hawking-Page transition is achieved for this model. Further, we observe that there exist stable black 
holes with small radii and that these regions are enlarged when choosing small values of the string cloud 
density and of the Gauss-Bonnet parameter. Besides, the rate of evaporation for these black holes are 
studied, determining whether the evaporation time is finite or not. Then, we concentrate on the dynam- 
ical stability of the system, studying the effective potential for s-waves propagating on the string cloud 
background. 



1 Introduction 

Black holes are one of the most exciting and fas- 
cinating objects on gravity physics, mostly be- 
cause their thermodynamical properties have a 
deep similarity with statistical mechanics J24), |7J, 
151 . In the last decades there has been a renewal 
of interest in higher dimensional spacetimes asso- 
ciated with this kind of objects. Moreover, a lot 
of effort to understand the quadratic corrections to 
the Einstein-Hilbert action has been put. 

At the present time, it is well-known that for 
a spacetime with dimension D > 5 the Einstein- 
Hilbert action of gravity admits quadratic correc- 
tions constructed from coordinate-invariant ten- 
sors which scale as fourth derivatives of the met- 
ric. In particular, when D = 5, the most general 
theory leading to second order equations for the 
metric is the so-called Einstein-Gauss-Bonnet (EGB) 
theory or Lovelock theory up to second order. This 
class of model for higher dimensional gravity has 
been widely studied, because it can be obtained in 
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the low energy limit of string theory (see [43 1). For 
a D-dimensional spacetime with D < 5, the Gauss- 
Bonnet (GB) term does not give any contribution to 
the dynamical equations. 

According to the recent theoretical develop- 
ments, a scenario in which the fundamental build- 
ing blocks of nature are extended objects instead of 
point objects should be considered quite seriously. 
In particular, one-dimensional objects (strings) are 
the most popular candidates. The study of the 
gravitational effects of matter in the form of clouds 
of both cosmic and fundamental strings has then 
deserved considerable attention (see |6|, [27J, |38|). 
Many authors have found exact black hole solu- 
tions in Einstein-Gauss-Bonnet models with dif- 
ferent kind of sources, for instance, with Maxwell 
and Born-Infeld energy-momentum tensor, matter, 
a radiating source and Yang-Mills fields (see j20] , 
1211 . 1291 , BD, (30), 1231). Moreover, black holes 
with a five dimensional Gauss-Bonnet term in the 
bulk and an induced three-dimensional curvature 
term on the brane were also found 1151 , 1161 . 

Also, some authors have put great attention to 
the thermodynamics of these solutions, as a way 
to characterize the black holes (see fll], ESL ITTBI , 
1171 , 1191 , |34l ). In particular, R. Cai has studied 
the phase structure of the topological black holes 
in the EGB theory with negative cosmological con- 
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stant (see |9)). More recently, in HU the authors 
have investigated the entropy of black holes in GB 
and Lovelock gravity using the Noether charge ap- 
proach. In addition, in 1321 the authors have ex- 
plored all the thermodynamics quantities for the 
EGB model with AdS backgrounds. They have 
concentrated on the role of the sign of the Gauss- 
Bonnet coupling constant a. 

We are interested in the construction of black 
hole solutions in a five dimensional spacetime 
with an energy-momentum tensor coming from a 
string cloud model. In the present work our start- 
ing point is to solve the full non-linear Einstein- 
Gauss-Bonnet equations for a static and spheri- 
cally symmetric metric. We study not only the 
structure of the geometry but also the thermody- 
namical stability of the system. Later, we shall 
concentrate on the dynamical stability, i.e. we will 
study the effective potential arising from the scalar 
waves perturbation satisfying either the massive 
or non-massive Klein-Gordon equation. 

The paper is organized as follows. In Section 
|2] we recall the basic facts of the Einstein-Gauss- 
Bonnet theory and of the string cloud model stud- 
ied by P. Letelier. In the next section we present the 
exact solution for a EGB theory on a static spher- 
ically symmetric 5-dimensional spacetime with a 
source provided by a string of clouds. As a 
consistency check, we reobtain the well-known 
Schwarzschild 5-dimensional spacetime from our 
solution as a particular case. Moreover, we also 
study the possible event horizons and singulari- 
ties. In Section |4] we focus on the thermodynam- 
ical properties of the aforementioned spacetime. 
More concretely, we obtain the Hawking temper- 
ature, the Helmholtz free energy, the heat capacity 
and the entropy. In particular, the entropy of the 
black hole was computed using the Wald prescrip- 
tion. The sign of the heat capacity allows us to de- 
termine the local stability of the system, whereas 
the sign of the free energy tells us if the system is 
globally preferred or not. 

In Subsection |4.3[ following th techniques of the 
tunneling method considered in |4], we compute 
the quantum corrections to the entropy, and we 
prove that it coincides with the standard expres- 
sions (cf. eq. (27) in (42]) under fairly more general 
assumptions. In the next subsection we explicitly 
determine when a black hole evaporates in a finite 
time, generalizing previous considerations in the 
literature. 

Finally, in Section [5] we consider a scalar field 
propagating in the bulk given by the spacetime 
and we determine the effective potential in this 
case and its third order approximation. 



2 Generalities 

2.1 The Einstein-Gauss-Bonnet theory 

We first recall the general facts about the Einstein- 
Gauss-Bonnet model on five dimensions. It is 
a generalization of Einstein Relativity, since its 
action also involves quadratic corrections con- 
structed from coordinate-invariant tensors which 
scale as fourth derivatives of the metric (see l28l ). 
The action for this theory coupled on a 5-di- 
mensional Lorentzian manifold (M, g) with matter 
fields is as follows 

S=l f v / ^g(^-2A + ai?2 ;B )d 5 x + 5 mat , 

where S ma t denotes the action associated with 
the matter fields, a is the so-called Gauss-Bonnet 
coupling constant, which we consider to be non- 
negative, g is the determinant of the metric g and 

c>2 d T?f J - u P a l E>2 a r> pfif 

-"-GB — K-livpoK- r + -ft — iH^Hr . 

Hereafter we will adopt the signature 
( — ,+,+,+,+) for the metric and we shall follow 
all the main conventions of General Relativity 
(GR) explained in 11371 . 

The corresponding Euler-Lagrange equations 
(also called EGB equations) for the metric are 

Qp. v = + Q MV + Q MV = T^, (2.1) 
where Tju„ is the energy-momentum tensor and 

+ 47? E pa — 2 7? 7? pcr7 M 

We point out that, as in the plain General Relativ- 
ity, the divergence of the EGB tensor vanishes. 

As a remarkable property of this theory we may 
say that the resulting equations of motion have no 
more than second order derivatives of the metric 
and the theory has been shown to be free of ghosts 
when it is expanded around the flat space, avoid- 
ing any problem concerning unitarity (see |8|). 
Further, it has been argued that the Gauss-Bonnet 
term appears as the leading correction to the effec- 
tive low-energy action of the heterotic string the- 
ory (see ED, 153] )- 

2.2 The string cloud model 

In this subsection we shall briefly review the the- 
ory of clouds of strings (see 11271 ). This cloud of 
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strings will be the source of the subsequent fam- 
ilies of spacetimes satisfying the Einstein-Gauss- 
Bonnet model. 

The Nambu-Goto action of a string evolving in 
the base manifold (M, g) is 



S 



NG 



-hd\°d\ 



where p is a positive constant related to the string 
tension, (A , A 1 ) is a parametrization of the world 
sheet E described by the string with induced met- 
ric given by 

dx^ dx v 

hab = 9,M/ Wd>fi 

and h is the determinant of h. 

Alternatively, we may describe the world sheet 
by a bivector E € T(A 2 TA/) of the form 

, dx^ dx v 



ab 



d\ a d\ b ' 



(2.2) 



where e ah denotes the two-dimensional Levi- 
Civita tensor given bye 01 = — e 10 = 1. The former 
satisfies the following identities 



EA1 [a s /3 7 ] = Qj y^E^E^ = 0, 



(2.3) 



where the square brackets indicate antisym- 
metrization in the enclosed indices. By the Frobe- 
nius' theorem, any bivector satisfying the previous 
identities leads to a parametrized surface such that 
(IZ21 holds. 

We would like to point out the following iden- 
tity 

E' i<T E CTr E T!y = hE 1 ^, (2.4) 

which will be used in the sequel and which follows 
directly from the definition of h and E. 

In this description, the Lagrangian density be- 
comes 

1 \ 1/2 

2 H ' 

Taking into account that = 2dC/dg^, we 

may obtain that the energy-momentum tensor for 
one string is given by T^ v = pE^E^/V^h. 
Hence, we shall consider the following energy- 
momentum tensor 



C=p(- -E^E, 



E^E/ 
V^h 



for a string cloud with density p. 

By making use of the Bianchi identity and the 
fact that the divergence of the Einstein tensor van- 
ishes we see that (-h)- 1 V M (pE^ <T )E ff !y E„ T = 0. 
Contracting the previous identity with E T „ and us- 
ing equation (2^, we obtain V^pE^E/ = 0. 



Employing a system of coordinates adapted to the 
parametrization of the surface, we see that in fact 



3 Exact solutions 



o. 



(2.5) 



From now on, we shall restrict ourselves to the 
case of static spherically symmetric 5-dimensional 
EGB spacetimes (M,g), so we may suppose that 
(locally) there exists a system of coordinates 
(t, r, <fi, 9, £) such that the metric g is given by 



g^dx* 1 ®dx u = -G{r)dt®dt- 



G{r) 



dr(&dr + r us, 



where G is a real-valued non-vanishing C°° func- 
tion of r and u is the volume form of S 3 for the 
remaining coordinates, namely, 

u = d£ <g> d£ + sin 2 (£)(d0 ® d6 + sin 2 (9)d<j> ® dtji). 

We shall also assume that the string cloud is 
spherically symmetric, so the only non-vanishing 
component may be E tr = — E rt , because E is a 
bivector. Therefore, the non-zero components of 
bivector are given by 



E^ 



E i(r)(*V\- «Vi). 



This source corresponds physically to a collec- 
tion of strings which are extended along the radial 
direction, distributed uniformly over the sphere. 
That is, the worldsheet of each string covers the 
t-r plane, which is the simplest case. 



3.1 Computation of I) 



fJ,V 



Since the only possible non-vanishing component 
of the bivector E is E tr = — E rt , we have 



^=E tr E * 



Taking into account that h = ±2E* r E tr = -(E tr ) 2 , 
we conclude that T\ = -p|E ir |, T r r = -p|E tr |, 
and the other components vanish. Furthermore, 
the identity J2.5II yields d r (r 3 T\) = 0. As a conse- 
quence, 

r r r = n = -^Q (3-D 

for some real constant a, that we will take to be 
non-negative for physical reasons. 



1 In our convention, we easily see that that the value of the 
tensor-energy momentum at a timelike vector (e.g. Ttt outside 
the black hole or T rr inside it) is always positive (cf. |37j) 
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We shall now discuss the energy condi- 
tions for this source. For the the energy- 
momentum tensor written in diagonal form as 
T^ v = diag(— cr,p r , 0, 0, 0), the weak energy condi- 
tion (WEC) means that a > 0, and p r + a > 0, 
whereas the dominant energy condition (DEC) says 
that a > 0, and \p r \ < a. Moreover, the strong 
energy condition (SEC) states that a + p r > and 
3p r + a > (cf. (36)). The physical meaning of a 
and p r is the energy density and the radial pres- 
sure, respectively. It is important to mention that 
the WEC implies the so called null energy condition 
(NEC) which says that a + p r > (cf. |3%1 ). 

In our model we obtain that a = a/r 3 and p r = 
—a/r 3 . We see thus that the NEC is always satis- 
fied. Also, the WEC and the DEC are satisfied for 
a > 0, otherwise these conditions are violated. On 
the other hand, we find that the SEC is violated for 
a > 0, whereas for the opposite case the SEC holds. 



3.2 Computation of the metric g ab 

The diagonal EGB equations il.l} for the energy 
momentum tensor given in the previous subsec- 
tion are 



Q\ = Q r r = -^, 9\ = 9 e e = 9% = Q. (3.2) 



After a rather long manipulation we obtain that 
the components of the EGB tensor are given by 

Q\ = —j (2r 3 A + 6rG - 6r + 3r 2 G") 



+ $raG'(l-G), 



+ 2aG"(l-G), 



G + 2rG' - 2aG' 2 



and the remaining EGB equations are trivial. We 
point out that the second family of differential 
equations | |3.2| | are a consequence of the first family 
in i3.2i . In consequence, G(r) satisfies the identity 



G - 2r 2 G + — 
Qa 



r 4 A 

3r 2 - 2ar + b — 



0, 



for some real constant b. Finally, there are two 
branches for the solution 



G(r) = l + -(1±VA 



(3.3) 



where 



16aa 



x , 8a(2a-36) 4 A 
A = 1+ Or* + 3 QAl 3.3 



Taking into account that a realistic physical solu- 
tion must become the classical Schwarzschild solu- 
tion of General Relativity in five dimensions with- 
out sources (i.e. for A = and a = 0) when con- 
sidering the limit a — > 0, we conclude that the so- 
lution 03.31 1 with plus sign has no physical interest, 
since in that case the mentioned limit does not ex- 
is^ Moreover, from the comparison of the expan- 
sion up to first order terms of the square root of the 
solution (13.31 and the classical Schwarzschild solu- 
tion of general relativity (with a = and A = 0) we 
conclude that 



Warn 



4 . 
-OA- 



16aa 
3r 3 



(3.4) 

In general, the solution may have two possible 
singularities: the usual singularity at r = and 
also the so-called branch singularity at r& > which 
is defined as follows. Since G(r) must be real- 
valued, we should only consider the domain val- 
ues r such that the radicand in d3.4b is non-positive. 
If 1 + 4a A/ 3 > this radicand is always non- 
negative, but for 1 + 4a A/3 < there exist values 
of r such that it is negative. If 1 +4aA/3 < holds, 
we define ?*& as the minimum positive real root of 
the quartic polynomial 

rt (1 + t aA ) + 16a(ar 6 + m) = 0, (3.5) 

if it has three different positive roots, or as the 
maximum positive root, otherwise. Hence, the 
metric is a priori only defined for r E (0,rb). So, 
as we said, if 1 + 4a A /3 < 0, there is a finite radius 
a singularity and the allowed domain of the radial 
coordinate is then no more < r < oo (cf. 1391 ). 

On the other hand, if r —> 0, G(r) approaches 
to 1 — y/m/a, which can be interpreted as the fact 
that the EGB term removes the metric singularities 
at the origin (cf. BT1 ). 

In a similar way, for the vacumm state (i.e. m = 
0, a = 0) for (13.3b the metric takes the form 



G(r) 



1 



1 



1 



-a A 



2a V V 3 

In order that G(r) is real-valued, we must have 
that aA > -3/4. Hence, if A > this condition 
is fulfilled since a > 0, whereas for a negative cos- 
mological constant A, the GB coupling constant a 
must be in the interval [0, — 3/(4A)]. 

If 1+4qA/3 < 0, then the expansion of G(r) near 
n is 



G(r) 



1 



4a J a 



2 A detailed analysis of the positive branch with a = 0, and 
A = was made in 1 31 1 
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Let us now suppose that 1 + 4a A /3 > 0. We may 
analyze the asymptotic behaviour of the previous 
expression in the limit r — > +oo, which is of the 
form 



G(r) 



A 



3(1 + J1+S4 



This leads to the well-known fact that the Gauss- 
Bonnet term provides a correction to the original 
cosmological constant A (cf. HO)- 



3.3 Event horizons of the black hole 

In what follows we are going to briefly discuss the 
general properties of the static solutions found in 
the last section. We have there studied the condi- 
tions for the function G(r) to be real-valued, but in 
order that it defines a metric it must be non-zero. 

By definition an event horizon is a hypersurface 
of the form r = in (M, g) where is determined 
as the maximum positive root of G(r\) = 0. There- 
fore, using (13.41 , rh is given by 



A 4 1 2 a 

y^h - + 3^ + (m - a) 



0. 



(3.6) 



From the considerations of the previous subsec- 
tion, in the case 1 + 4a A/3 < 0, there may be an- 
other constraint on the existence of an event hori- 
zon, indicating that the event horizon should also 
satisfy the inequality < r & . We shall assume that 
this is always the case if we suppose that an event 
horizon exists and the inequality 1 + 4a A/3 < 
holds. 

If A = 0, then equation (13.61 becomes quadratic, 
and its solutions are 



Hi 



= a ± \f a 2 + 2(m — a), 



(3.7) 



where a = a/ 3. In this case the black hole solution 
has an inner horizon (smaller positive root) and an 
event horizon (largest positive root or outer hori- 
zon). This solution exhibits a similar causal struc- 
ture of charged Reissner-Nordstrom (RN) black hole 
geometry in 5 dimensions. As a consequence of 
eq. (13.71 1, the horizons for A = are well-defined if 
and only if 



2(m 



> 0. 



The special case where the equality holds is re- 
ferred to as an extreme black hole with the event 
horizon given by r cxt = a. If a 2 + 2 (m — a) < 0, the 
string cloud solution describe a naked singularity 
at r = 0. 



Now, we establish what the leading terms of the 
most relevant quadratic invariants near the physi- 
cal singularity (r = 0) are: 



bed 



RabcdR"" 
K a bH 



12m 4a 



12m 



a r 
6a 

ar 3 



0(r- 2 ), 
0(r- 2 ), 



It turns out that the mass term dominates over 
the string cloud contribution near the singular- 
ity. However, the last term is smoother than 
the one for the 5-dimensional Schwarzschild met- 
ric with negative mass, that possesses a timelike 
naked singularity at the origin with Kretschmann 
scalar R a bcdR abcd diverging as r" 8 (cf. [21 1). On 
the other hand, the Kretschmann scalar for the 
charged Gauss-Bonnet black hole has a leading 
term of the form r~ 6 but it also has similar terms 
as in our case. 



4 Thermodynamics of the black 
hole 

In this section we shall discuss and reckon the 
main thermodynamical properties of the string 
black holes solution within the EGB framework. 
For the entire section a will be treated as an exter- 
nal parameter in relation with the thermodynami- 
cal considerations. 

Using the standard definition of surface gravity 
AC = l/2|G'(rh)| for a spherically symmetric static 
metric with a static Killing horizon, it follows that 



AC = 



I'-h 



A 31 
3 Mil 



rfj + 4a 



Using the identity Th = AC /2n, we see that 



T H = 



A r 3| 
3 r hl 



2ir(r 2 + 4a) 



(4.1) 



(4.2) 



Since it will useful later, we remark that in the for- 
mula of the Hawking temperature we have omit- 
ted the Planck constant h in the numerator by our 
choice of units. 

Let us consider the case with vanishing cosmo- 
logical constant. In the limit a — > and a — > 0, 
we recover the classical result for Schwarzschild 
spacetime saying that the Hawking temperature 
diverges when the event horizon shrinks to zero. 
On the contrary, for vanishing a but non-zero EGB 
coupling constant a, the Hawking temperature re- 
mains finite when performing the limit — > 
(cf. |9) )■ Besides, for a/0 and a^O, there is a par- 
ticular radius of the event horizon for which the 
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Hawking temperature vanishes, so the black hole 
does not radiate energy. As a general behaviour, 
the Hawking temperature exhibits a peak, which 
decreases and moves to the right when the string 
cloud parameter a grows (see Fig. [TJl. 



0.0(1 









a=),A=-o.l 




.a=0.001 
//a=0.01 ^^^^ 
// _,-<S=0. 1 
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Figure 1: The Hawking temperature as a function 
of the horizon radius for A = —0.1, a = 1 and 
different values of a. 



4.1 Global stability 

It is a well-known fact that asymptotically 
AdS Schwarzschild black holes are thermally fa- 
vored when their temperature is sufficiently high, 
whereas the pure AdS background is thermally 
preferred for low temperatures. 

As the temperature decreases, there is a first or- 
der transition such that the black hole spacetime 
leads to a pure AdS geometry. This effect is known 
as the Hawking-Page transition (HPt) (see [25]). The 
phase transition has been widely studied for oth- 
ers asymptotically AdS black holes in the context 
of EGB gravity (cf. E), El, (TBI, EL El, EH)- 

Now, we are going to study the thermodynamic 
phase structure for EGB gravity with a source 
given by a string cloud. In order to do so, we 
first need to express the black hole mass param- 
eter m as a function of the event horizon radius rv 
which can be done using J3.6I >- The mass parame- 
ter m is related to the ADM mass M of the black 
hole, which takes the form M = 'iVt^m/ (87r) for a 
5-dimensional spacetime, where VI3 is the volume 
of the unit sphere in M 3 . 

Our next step is to compute the entropy in the 
case A < 0. So, from the first law of thermody- 
namics dAI = TdS and the fact that A < 0, we can 
obtain the entropy of the black hole 

s= r+dM dr+ 
Jo dr + T(r+) 1 

where r + denotes the radius of the largest horizon 
event (see El)- We have imposed the physical as- 
sumption that the entropy vanishes when the hori- 



zon of the black hole shrinks to zero as explained 
by Cai (see [9 1). The expression is given by (cf. |3T1 ) 

S = tt(^± +Aar + y (4.3) 

As it is well-known, the entropy for the EGB model 
does not satisfy the Bekenstein-Hawking law S = 
A/4, since S = A(l + 12a/r\)/A, where A = 
47T7+/3 is the area of the (spherical) event hori- 
zon of radius r + (cf. |9|). Interestingly, the energy- 
momentum tensor of the string cloud does not 
lead to any correction in the entropy d4.3l) . 

In the previous computation we have asssumed 
that the ADM mass M is the energy of the black 
hole. In order to prove so, we might proceed as fol- 
lows: first, compute the entropy function S\y us- 
ing the Wald formulation and then, from the first 
principle of thermodynamics dSw = MyfdTs, de- 
duce the (real) mass of the black hole A% . We re- 
call that the Wald formulation gives an expression 
of the entropy for any A^-dimensional spacetime 
with a diffeomorphism invariant Lagrangian such 
that it admits stationary black hole solutions with 
a bifurcate Killing horizon (with bifurcation sur- 
face £), and that the canonical mass and angular 
momentum of the solutions are well defined at in- 
finity (see 1401 ). We shall apply this formulation 
to the spacetime we are considering, for which the 
entropy is given by 

S w = -2tt / — e ab e cd Vhdn 3 , (4.4) 

j£ 0-tlabcd 

where C = £egb + £mat is the sum of the La- 
grangians corresponding to the EGB model and 
the matter source, respectively, \fhdri3 is the vol- 
ume element induced on £ by the metric of the 
spacetime, and e a b is the binormal tensor of the 
spatial section £ defined by the horizon given by 
r = rh and t = etc, satisfying the normalization 
condition e a be ab = —2. Giving two null (local) vec- 
tor fields £ and v normal to S and satisfying that 
(?(£, v) = 1, the binormal tensor may be written as 

£ab = ia^b — Va£,b- 

We point out that in our case SC m m/SR a bcd=0, 
due to the fact that the string cloud matter is min- 
imally coupled to gravity. Therefore, the matter 
term in (|4.4b does not contribute to the entropy 
Sw- It is not difficult to see that for the variation 
of the gravity action with respect to the Riemann 
tensor, regarding R abc d as formally independent to 
the metric g a b, coincides with the entropy given in 
d4.31 >- A detailed derivation of this result can be 
found in Q or OS. 

Moreover, using the Wald's entropy and the first 
law of thermodynamic dAfyj = TndSw we get, up 
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to a constant Co, the mass of the black holes as 

2M w (r + ) = m(r+) + C . (4.5) 

Taking into account that we should obtain the 
Schwarzschild spacetime in the GR limit stud- 
ied in Subsection 13.21 we conclude that the pre- 
vious mass 2Mw is equal to the mass of the 
Schwarzschild black hole, which further implies 
that Co vanishes. Hence, Mw coincides with the 
ADM mass M. 

The free energy F = M — T^S of the black hole 
is given by 

_ + 2(1 + QaK)r\ - 8ar 3 + - 12ar 2 + + 48a 2 

2A(r\ + 4a) _ ' 

(4.6) 

For large values of r+, F diverges to sgn(A)oo. On 
the other hand, for small values of r+, F goes to a. 

In the limit a — >• 0, Il4.6b gives the same expres- 
sion found in |9j. For A < we obtain that for 
small radius of event horizons the free energy is 
positive, so the black hole is globally unstable (in the 
thermodynamical sense), whereas for large values 
of r + it is negative {i.e. globally stable black hole) 
(cf. |3l). This is thus the Hawking-Page transition 
mentioned above. Besides, using the eq. 114.61 , we 
can explicitly obtain the critical temperature where 
the free energy is zero (F(T C ) = 0). Then, the criti- 
cal temperature in term of the parameters A, a and 
the outer horizon r + is given by 

24tt(^- +4ar + ) 

The latter critical temperature essentially indicates 
the point where the HPt occurs. More precisely, for 
T > T c the black hole solution is thermally glob- 
ally preferred with respect to the reference back- 
ground while for T < T c the reference background 
solution is globally favored (for a detailed discus- 
sion on this point with a = and with Maxwell 
corrections, see |3l). 

3 — , , — , , — , , — , , , — , , — , , — , , — 

A=-0.< a=0,a=O. 

2 \ : 

1 \ : 

_^ -0=0.5,3=0.— «. \ - 

o — -^^"rrrrrrrr — v 

F , \ a=0.5, a=J>^ \ 5 

— 1 \ s \ 

\ \ 

-2 \ \ \ I 

a=0.5,a=2>*. \ \ 

-3 \ \ \ : 

_ 4 t — , , — , , — , , — , , , — , , — , , — V V 

2 4 6 8 

r+ 

Figure 2: The free energy as a function of r + for a 
fixed value of A and several values of a and a. 



In Fig. |2]we indicate how the HPt is achieved by 
varying the value of the string cloud parameter a: 
if the free energy only admits one zero the system 
passes from an unstable phase to a stable regime. 
However, from a more careful analysis we see that 
the free energy may have three roots, so the system 
exhibits two stable phases: one for the small radius 
and the other for the large horizon radius. Besides, 
we obtain that the first stability zone is enlarged 
when increasing the parameter a over the range 
[0.5,1.3]. 

In Fig. [3] we fix A = —0.1 and a = 1. We ex- 
amine F(r + ) when the Gauss-Bonnet parameter 
varies over the range [0.001, 0.43]. It turns out that 
the function F{r+) has the same phase structure, 
U-S-U-S (U denotes an unstable zone and S a stable 
one), then our model seems not only to present the 
standard HPt but also to offer new phases for the 
global stability of the EGB gravity. Moreover, the 
first stability region corresponding to small hori- 
zon radii is amplified when increasing a. 

Here we have only dealt with the case a > be- 
cause our model is more related to string theory. 
However, some authors have considered the EGB 
model for a < and its implications in the ther- 
modynamics of black hole solutions (see 1 32 1). 



a=1. A=-0.7 






n=0.001,/ / \ 4 








/ / / ,»«=a3>^-- N \ \ \ 
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Figure 3: The free energy F as a function of r + for 
a fixed value of a = 1 and A = —0.1. 



4.2 Local stability 

Despite the fact that we have studied the condition 
for the global stability of the solution we are also 
going to examine the local thermodynamics struc- 
ture by computing the specific heat. The reason is 
that even when a black hole configuration is found 
to be globally stable, it could also be locally unsta- 
ble to some globally favored configuration (cf. ||3l 
and references therein). 
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The heat capacity C = dM j&T turns out to be 

c _ 7r(rg + 4 Q ) 2 |r h -g- A r 3| 

(1 - Ar 2 )(r 2 + 4a) - 2r h (r h -a- f r$) ' 

(4.8) 

By simplicity, we shall only analyze the sign of C 
for A = 0. For a — > we recover the specific heat 
formula reported in [9[. 

It is a well-known fact that the local thermody- 
namic stability of the system is related to the sign 
of the heat capacity. When rh > a, we get that the 
black hole is locally stable to thermal fluctuations 
{i.e. C > 0) iff m < a, whereas the heat capacity 
is negative (so the black hole is locally unstable) iff 
m > a. Notice that the point where the heat capac- 
ity vanishes is a = m, corresponding to rh = 2d. 
Furthermore, in this case, the critical temperature 
is 



M 87r(a 2 +a)' 

In the opposite case rh < a, the thermodynami- 
cal behaviour is interchanged, i.e. the solution is 
stable (resp. unstable) iff m > a (resp. m < a). 
The jump of the heat capacity at the point rh = 2d 
indicates a second order phase transition (see Fig. 
IDl. Notice that in the GR limit (a -> 0), the critical 
temperature remains finite. 



CO 





su=0.5, 3=0.5 
_,,-'*' s a=0.9, a=0.02 


\ S o=0.5 


a=Q$ 

\ ✓ a=0, a=0 




! a=0.5, a=0.5 s \ 


=0, a=0s\ 


A =-0.1 



2 4 6 8 



Figure 4: Heap capacity C as a function of r + for 
A = -0.1 with (a, a) = (0, 0), (a, a) = (0.5, 0.5) or 
(a, a) = (0.9,0.02). 



4.3 Entropy: Quantum corrections 

In this subsection we shall analyze the quantum 
corrections to the classical entropy computed in 
J4.3I > for A = 0. Since we shall consider formal se- 
ries on the Planck constant h, we shall restore it in 
the previous formulas for the entropy. 

We shall denote 5bh = A/ (Ah) the entropy 
given by the Bekenstein-Hawking law and S*gb the 
entropy computed in the previous section (with 



the h) 

A ( 12a \ 

In the semiclassical approximation for General 
Relativity, the black hole entropy satisfies the 
Bekenstein-Hawking area law, whereas when the 
full quantum effect is considered, this area law 
should be corrected. Further, the corrected en- 
tropy 5 q takes the form 

S q = S + 9\n(S) + (4.9) 

where S is the uncorrected semiclassical entropy 
of the black hole and 9 is a dimensionless constant. 

Following the considerations on the tunneling 
method of [4], the corrected Hawking temperature 
T q of a black hole is given by 

i>l 

where Th is the standard semiclassical Hawking 
temperature of the black hole and for some (di- 
mensional) constants 7; . 

Whereas in the first three paragraphs of this sec- 
tion we have considered a and a as external pa- 
rameters, we shall now allow them to vary within 
the model. This in turn implies that they should 
have a corresponding work-term in the expression 
of the first of thermodynamics, i.e. 

dM = T H dS G B + W a da + W a da, 

where W a and W a are the work-terms of a and a, 
respectively. At this point, one may argue why 
these entropy and temperature should coincide 
with the one we computed in the previous subsec- 
tions. We leave the proof of this simple fact to the 
Appendix. 

Now, we see that, since the work-terms of a and 
a are not supposed to change by quantum effects, 
the semiclassical approximation leads to 

T H dS* G B = T q dS qi (4.10) 

which means that in fact we need not compute the 
work-terms coming from the parameters a and a. 

As noted in [42] , the constants 7, should be 
determined using dimensional considerations. In 
this case, since the following dimensions coincide 
[a] = [a] 2 = [r] 2 , we propose 

Cti 

"~y ' = r 

/v^ i 71 -J 3-Q1+J2) y ' 

( Yj a -- r + ana 5 

where a, and a-, are dimensionless constants and 
the sum is indexed by j = (ji, j'2) € M 2 and of 
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finite support. This is a generalization of the eq. 
(27) proposed in l42l . 

Hence, the corrected Hawking temperature has 
the form 

I v 1 

Ea.j K 


i>l Va'rfa^a 2 

V ~ { G) 3 + J J 

(4.11) 

From equation d4.106 , the differential of the cor- 
rected entropy is given by 



dS a 



((7rr + + 4ira)dr + + 47rr_|_da), 



where we have used eq. [ 14.31 1. Since the entropy S q 
is a state function, dS q should be an exact differen- 
tial, so 



d ( T n, 2 „ \ 



do 

d /T H 



d /T H 



di 



(f^ +4 -))=°'l(| 4 -)=° 



dr. 



From the last two equations we see that Tn/T q 
should not depend on a, so 



a,- ft 4 



V 



(ie») 



(4.12) 



The remaining equation implies that 



(J"6R) 



i-2 r+ a 2 



Hence, for all j e M we have that 12 (j — l)a* = 
(5 — j)a l j-2- Since the support of the sum was fi- 
nite, we see that a* = for all j < 1 and j > 5. 
Moreover, by the form of the previous expression 
we conclude that a* vanishes it j ^ 2Z + 1 and if 
j = 5. The only non- trivial relation we obtain is 
12a| = a\, which means that 



7r(ri + ^ra)* ' 



for some new dimensionless constant ati (cf. Il42l , 
eq. (27) and eq. (32)). So, 



i>l GB 

which in turn implies that (using eq. 1 14.101 )) 



5q — SqE 



log(5 GB ) + ^ 



(1-0Sgb 



const. 



This expression coincides with the one found in 
eq. (27) in [42], which seems reasonable since the 
entropy does not depend on the extra parameter 
a. However, we would like to stress that we have 
adopted a more general assumption than 1421 for 
the expression of the constants in the semiclassical 
expansion of the Hawking temperature. 

4.4 Evaporation rate of the black hole 

We shall discuss the evaporation of the black hole 
with A = 0. Taking into account the Stefan- 
Boltzmann law for the black body in 5 dimensions 
we see that the evaporation rate is 



dm 
~dt 



ex 



-Thr 



(4.13) 



First, we see that we have a well-defined horizon 
radius if and only the mass parameter fulfills the 
requirement 



> sup.{0,a - a 2 /2}. 



(4.14) 



This in turn implies that, if a — a 2 /2 > 0, then the 
black hole cannot decay (this generalizes previous 
considerations in |2]). 

Otherwise, if a — a 2 / 2 < 0, the black hole could 
decay or not. More concretely, if a — a 2 / 2 = then 
the evaporation time of the black hole is infinite, 
whereas if a — a 2 j 2 < it evaporates in finite time 
for a > and in infinite time for d = 0. This can be 
proved as follows. 

Since dm/dt = (r + — a)dr + /dt, we may assume 
that, as m varies from the initial mass parameter 
mo at time to to at time t e , r + goes from r + .o to 
r +-e = a + \Ja 2 — 2a. Then, if we denote by At = 
t e — to the interval of time for that variation 



At 



(£)"V 



oc 



r +,o 



a) 4 r 3 



dr. 



If a — a 2 j 2 = 0, then one limit of integration is 
r + . = a and the integral diverges since the inte- 
grand goes as (r + — a)~ 4 . Hence the evaporation 
time is infinite. 

On the other hand, if a—d 2 /2 < 0, and assuming 
that a > 0, we see that the integrand is bounded 
and continuous within the integration interval, so 
the black holes evaporates in finite time. In the 
case a = 0, we easily see that the integral diverges. 

5 Scalar waves propagating on 
string cloud background 

In this section we shall consider the case of a scalar 
field $ with mass ii on a 5-dimensional spacetime 
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with a string cloud background. We follow mainly 



The scalar field $ satisfies the Klein-Gordon equa- 
tion 

-^d a (^g ab d b <5>) - ,i 2 $ = 0. 



We shall assume that, in local coordinates, 
®(r,t,(j>,d,£) = *(r,6',^)e- ,; ( tJt - Zl «, where h £ Z, 
and further that *(r,6>,£) = R(r)Y(9)Z(£). Since 
the solution in the angular part is the same as in 
the classical case, we shall only analyze the radial 
part. In this case, there exist a constant I3 E No 
such that the radial Klein-Gordon equation can be 
rewritten as 



d r (r 3 Gd r R) fu 



rR 



(5.1) 



Now we choose R(r) = r 3 ^ 2 R(r) and the so 
called tortoise coordinate = r*(r) given by 
dr*/dr = G(r)" 1 . We also set R*(r*) = R(r). 
Therefore, if we divide equation J5.1II by r -1 / 2 , we 
get 



d 2 R* 
dr 2 



+ L 2 -V(r)^R* 



0, 



where 



/'" 



Using the expression for G(r) given in ( |3.4b and 
the dimensionless variables f = r j ' \fm, A = Am, 
a = 2a/m and a = a/^/m, we obtain that G(f) = 

1 + fsf 1 ~ y( f ))> where 



y{r) 



1 



2d A 8aa 8a\ 1 / 2 



3f 3 



Let us now assume that A = 0. This implies that 
there are no terms of positive exponent in the ex- 
pansion in f of the potential. Its third order ap- 
proximation is 



V(f) ~ /i 2 



2a/i 2 



H5 + !-V)^-^. 



We see that the second order approximation co- 
incides with the one in [26 1 (when A = a = \i = 0) 
and that the parameter a does not provide any cor- 
rection up to third order. We remark that the pa- 
rameter a appears in the fourth order term of the 
expansion. 

Moreover, from the analysis of Fig. [5) we see 
that the potentials are real and positive outside 
the event horizon. Hence, following the argument 
given by S. Chandrasekhar in 1121 , the black hole is 
stable under scalar perturbations. Besides, when a 



maximum of V(r) exists, it seems to decrease as a 
increases. 




Figure 5: Potential V as a function of r + for a unit 
mass parameter m = 1 with fixed values of A = 0, 
a = 0.5 and fx = 1 and different values of a and Z3. 



6 Conclusions 

To summarize, we have found black holes in 
Einstein-Gauss-Bonnet gravity for a static and 
spherically symmetric 5-dimensional spacetime 
with an energy momentum given by a cloud of 
strings. We characterized the solution, calculat- 
ing the possible horizons, which could be at most 
two. Besides, we confirmed the singular structure 
of the spacetime at the origin by computing some 
quadratic invariants. 

We performed a detailed analysis of the 
more relevant thermodynamical aspects, focusing 
mainly on the global and local stability of the sys- 
tem. In particular, the entropy of the black hole 
was computed using the Wald prescription. We 
also found that the HPt can be realized within this 
context. Interestingly, the region for small stable 
black holes are considerably enlarged when choos- 
ing small values of (a, a). This has been already 
observed in the literature for other sources (cf. 

m, ma, my 

Regarding the local stability of the black holes, 
we obtained that the heat capacity can be negative 
or positive definite, depending on the relation be- 
tween in and a, which further tells us whether the 
black hole is stable or not. In particular, the point 
where the heat capacity is zero corresponds to the 
case 7"h = 2a and non-zero Hawking temperature. 
Notice that this critical temperature is zero for a 
vacuum solution a = (cf. |8|, |9|). 

Later on, we examined the quantum correction 
to the Hawking temperature and to the entropy, 
following the prescriptions of tunneling method 
used in [4], [42]. We have provided a proof of a 
lemma which tells us that the expression for the 
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GB entropy found in the case that a and a are 
treated as external parameters holds when these 
parameters are allowed to vary. 

After that, we estimated the evaporation rate 
for the black hole , obtaining that for 2a > a 2 
the black hole cannot decay for its mass parame- 
ter is bounded from below. However, for 2a < a 2 
the black hole evaporates in finite time if a > 0, 
whereas the evaporation time is infinite if a = 0. 
For the critical case corresponding to 2a = a 2 we 
showed that the evaporation time is infinite. 

Finally, we studied the propagation of scalar 
waves using this metric as a background. We ba- 
sically found that when the potential V(r) has a 
maximum, the black holes are stable. 

Acknowledgments 

We are grateful to the referee for his careful read- 
ing of the manuscript. This work has been par- 
tially supported by the projects UBACyT X212 and 
PICT-2007-02182. The first author is a research 
fellow of the Alexander von Humboldt Founda- 
tion. The second author is a CONICET fellow (Ar- 
gentina). 

Appendix 

In this appendix we shall prove the elementary re- 
sult stated in Lemma 16.11 We may rephrase it in 
more general terms (e.g. for Banach manifolds) but 
we prefer to stay concrete. We first need to state 
some notation. 

We fix p > 1. Let us suppose that two C p func- 
tions / : U C E"+ m — > E and g : U C E ,l+m -> 
E n " = M„(R) defined on an open set U are given, 
such that det(g) is non-vanishing. We shall denote 
an element of R n+m as (x, y), with x £ 1" and 
y £ M m . 

We are interested in a solution to the following 
problem: find C p -morphisms / : V C E n+m -> K 
and h : U C R n+TO -> R™ defined on an open set 
V such that the map H : U -> V given by (x, y) H> 
(h(x,y),y) is C p -isomorphism and the following 
identities hold in U 

%(h(x,y),y)=g(x,y), 
f(h(x,y),y) = f(x,y). 

A problem as before will be called a change of do- 
main problem. Notice that (/, h) are not indepen- 
dent since / = / o H^ 1 . We also see that 

|r(»,y) = sfoj/) -1 !^^*/)- (6-2) 



Furthermore, the fact that H is a C p -isomorphism 
implies that the previous expression is an invert- 
ible matrix. 

Lemma 6.1. Let us suppose that two C p functions f : 

U C E n+m -> E and g : U C W l+m -> M" 2 = 
M n (E) defined on an open set U, are given, such that 
det(df/dx(xo,y~o)) an d det(g(xo,yo)) do not vanish 
for a fixed point (x ,j/o) £ U. Then a local solution 
to the change of domain problem in a neighbourhood of 
(xo,y~o) exists. Moreover, any two (global) solutions 
to the change of domain problem (fi,hi) and (f 2 , /12) 
satisfy that hi — h 2 is a function ofy £ E m . 

Proof. Using the Existence and Uniqueness Theo- 
rem for ODE's in eq. 116.211 , we see that there exists a 
locally defined function h as needed. The Inverse 
Function Theorem tells that the map H given by 
(x,y) 1 — y (h(x,y),y) is a local C p -isomorphism at 
(xo, yo). We restrict the domain in order to assure 
this last property. Now we define / = / o H~ l . 

We need to check that the first equation in Il6.1b 
holds. However, from the chain rule and the In- 
verse Function Theorem, it is easy to see that the 
first equation in | |6.1| | is equivalent to d6.2| |, and so 
the local existence follows. 

Also, we see that if a solution to the previous 
problem exists, it is not unique. This is due to the 
following simple fact: if k : W' 1 -> E" is a C p - 
morphism, we may define h(x, y) = h(x, y) + k(y). 
It is trivially verified that the map H given by 
(x,y) 1 — y (h(x,y),y) is a C p -isomorphism, being 
the composition of H and the C p -isomorphism 
(x,y) h-> (x + k(y),y) (with inverse (x, y) i-> (x — 
k(y),y)). We may now define f(x,y) = f(x — 
k(y), y). It is trivially verified that (/, h) is also a 
solution to the change of domain problem. More- 
over, given any two solutions to the change of do- 
main problem, they are related as before, since 
equation | |6.1| | implies that 

dh _ ,9/ 

g={x,y) = {g(x,y)) ^(^W- 

The lemma is thus proved. □ 
We may apply the previous lemma to the case 
(n, m) = (1, 2), x = r + ,y = (a, a), / = M(r+, a, a), 
and g = Tn(r + , a, a). In this case, it it easy to see 
that the function / = S given by eq. d4 3b (and 
a certain h) provides a solution to the change of 
domain problem. Furthermore, by the previous 
lemma and the physical assumption that the en- 
tropy vanishes for r + = we conclude that the 
expression for the entropy given by 04. 3t is unique 
and it also holds when allowing the parameters a 
and a to vary. 
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